Frontiers in Molecular Biosciences

Background: Citrin deficiency, caused by biallelic pathogenic variants in SLC25A13, must be identified early to prevent serious complications such as hyperammonemia and liver failure. However, clinical diagnosis is often delayed due to its nonspecific presentation and limited sensitivity of amino acid-based newborn screening methods. Although genome-based evaluations are being investigated to address these issues, concerns about their cost, turnaround time, variant interpretation ability, and data handling highlight the need for a more practical yet reliable alternative. We investigated the feasibility of applying proteomic approach on dried blood spots (DBS), which are routinely used in newborn screening. Methods: We performed untargeted liquid chromatography-tandem mass spectrometry to analyze the proteome of DBS using a previously developed "non-targeted analysis of non-specifically DBS-absorbed proteins" (NANDA) workflow. SLC25A13 protein abundance was quantified in individuals with biallelic loss-of-function mutations, compound loss-of-function/missense mutations, and heterozygous carriers; this was also evaluated in healthy and diseased controls representing relevant differential diagnoses. To leverage proteomic information, we derived a multivariate proteomic signature using feature selection and evaluated its performance with leave-one-out cross-validation. Biological relevance was assessed by enrichment analysis, and complementary transcriptomics was performed using RNA sequencing. Results: A total of 7,474 proteins, including SLC25A13, were consistently detected in DBS. SLC25A13 was undetectable in individuals with biallelic loss-of-function mutations. However, individuals with compound loss-of-function/missense genotypes showed reduced but measurable SLC25A13 levels, comparable to those observed in heterozygous carriers. In contrast, a compact 15-protein signature accurately identified individuals with compound loss-of-function/missense genotypes (AUC, 0.99; sensitivity, 1.00; specificity, 0.95). The signature was enriched for Ca2+-response, and transcriptomics showed downregulation of genes related to multimodal ion channels in affected individuals compared to controls. Conclusions: DBS-based proteomic profiling may assist in the diagnosis of citrin deficiency through SLC25A13-quantification and a biologically plausible multivariate signature. More broadly, this strategy offers a promising new diagnostic layer for protein disorders, providing a proteomic readout in a clinically practical DBS format with potential utility for future diagnostic and screening applications.

Biomolecular condensates that form via liquid-liquid phase separation (LLPS) of, most prominently, intrinsically disordered proteins (IDPs) are ubiquitous in eukaryotic cells and responsible for regulating a plethora of biological functions. Amongst these, they contribute to regulating cell motility, either individually within an extracellular matrix or collectively within confluent epithelial tissue. In this computational study we focus on the latter with the aim of investigating whether the mutual exertion of mechanical forces during collective migration in an epithelium can principally trigger cytoplasmatic LLPS. Since present models for confluent epithelial motility have so far only considered cells that are devoid of phase separating (protein) solutes, we extend a common multiphase approach for 2D cell motility with a mixing contribution including any number of protein solutes. Our model considers the phase behavior in both intracellular and extracellular regions and determines to what extend the membrane is permeated by the solutes under the influence of mechanical and osmotic forces. Our initial calculations unlock a very rich behavior involving formation and dissolution of condensates during migration, as well as an impact of LLPS on the very nature of the motility itself, through feedback mechanisms which may bear biological relevance.

Coarse-grained structure-based models (CG-SBMs; or G[o] models) are simplified potential energy functions of biomolecules or biomolecular complexes that encode their structure. Molecular dynamics simulations of such SBMs have been successfully used to study long time-scale dynamics such as protein and RNA folding, and large conformational transitions of biomolecular complexes. SBMs have several advantages: (1) Their MD simulations are computationally inexpensive, making extensive sampling easily accessible to many researchers. (2) They are easy to modify and can be adapted for the specific biomolecular problem that needs to be investigated. However, the force-fields of SBMs are not usually included in commonly used biomolecular simulation packages resulting in a barrier to their use. Here, we present SuBMIT (Structure Based Models Input Toolkit; https://github.com/sglabncbs/submit), a toolkit for generating coarse-grained SBM input files for performing MD simulations with GROMACS and OpenMM/OpenSMOG. Simulations whose input files can be generated using the different flavors of CG-SBMs present in SuBMIT include the folding and conformational ensembles of proteins with intrinsically disordered regions, 3D-domain-swapping in proteins and the dynamics of RNA-protein assemblies (e.g., simple RNA viruses).

Surface factors that contribute to the virulence of Staphylococcus aureus have become therapeutic targets in the treatment of illness associated with this bacterium. Staphylococcal protein A (SpA) is a well-known contributor to S. aureus toxicity and virulence, although relatively little is known about protein A and how its biological function has evolved. SpA is displayed on the surface of the bacterium and contains 5 nearly identical helical ({approx} 60 aa) domains that bind antibodies with high affinity (Kd {approx} 10 nM). The folding free energies of only domains E and B have been determined. In this study we used intrinsic fluorescence detected denaturation to measure the folding thermodynamics of each domain in isolation and in the native multidomain context using a construct that includes the N-terminal half of the mature protein (SpA-N). We also constructed a series of proteins with 1 to 5 repeats of B domain, linked exactly as the five domains of WT SpA are linked. We used nearest neighbor thermodynamic models to explicitly demonstrate that the domains in B domain repeat proteins fold independently. We also showed that the domains in SpA-N fold independently by comparing the folding free energies of domains in isolation and in their multidomain context. Previous dynamic NMR experiments detected highly flexible linkers between domains in 5B, suggesting that the domains of SpA are structurally independent, which is likely responsible for the lack of thermodynamic coupling. Our results also showed a steep increase in domain stability from the N-to C-terminus in SpA-N, from 0.97 {+/-} 0.05 to 5.57 {+/-} 0.28 kcal/mol. We hypothesize that this stability gradient is related to efficient secretion of protein A.

The Na+/K+-ATPase (NKA) regulates ion balance in the kidney and influences cellular processes like proliferation and apoptosis through its signal transduction. The endogenous ligand 20-Hydroxyeicosatetraenoic acid (20-HETE) contributes to inflammation and fibrosis in chronic kidney disease (CKD) and inhibits NKA activity in renal tubules. However, the molecular mechanism of this interaction remains unclear. In this study, we used in-silico approach to investigate the potential interaction between 20-HETE and NKA. Various ligands, including known NKA ligands such as cardiotonic steroids (CTS), 20-HETE, and negative controls, were docked using rigid and Induced Fit Docking to predict the affinity of the ligands toward NKA. Binding free energy calculations with the Prime Molecular mechanics with generalized Born and surface area (Prime MM/GBSA) tools were used to confirm the involvement of key amino acids in ligand-receptor interactions. The docking analyses revealed that 20-HETE exhibited a binding affinity comparable to negative control, with some differences between rigid and induced fit docking. Binding free energy data highlighted key amino acids in the 20-HETE and NKA interaction. Interaction fingerprint and mutations such as Ala330Gly and Val329Ala significantly reduced binding free energy, while Thr804Ala showed a notable decrease, underscoring the potential importance of these amino acids in ligand stabilization. These findings provide computational evidence supporting potential direct interaction between 20-HETE and NKA and identify candidate residues for future experimental validation.

Mammalian single-stranded DNA binding protein 1 (SSB1) has been established as an essential component of genome stability in both human cells and mice. Moreover, SSB1 was recently implicated in cytoplasmic stress response by its involvement in Ras GTPase-activating protein-binding protein 1 (G3BP1)-containing cytoplasmic stress granules (SGs) upon various forms of stress. Here, we generated and analyzed human cellular knockout and rodent ischemia-reperfusion (I/R) models to define SSB1s roles in cytoplasmic stress response. Analysis of wild-type as well as SSB1 and G3BP1 knockout human retinal pigment epithelial (RPE-1) cells shows stress-specific incorporation of SSB1 into SGs and a negative regulator role for SSB1 in SG dynamics under sublethal stress conditions. We find that SSB1 knockout measurably increases cellular sensitivity to oxidative stress but does not alter cell proliferation following mild acute stress. Moreover, we detect SSB1 efflux from the nucleus upon stress that is dependent upon the presence of G3BP1 in a stress-specific manner. In addition, using mouse and rat models we observe significant upregulation and robust cytoplasmic granulation of SSB1 upon renal ischemia-reperfusion stress, establishing SSB1s involvement in complex organismal stress response in vivo. Together, our data demonstrate active involvement of SSB1 in cytoplasmic response to cellular stress and acute kidney injury, with implications for targeting stress response functions in cancerous versus non-cancerous contexts. HIGHLIGHTSO_LISSB1 is incorporated into cytoplasmic stress granules and negatively regulates stress granule assembly under sublethal stress conditions C_LIO_LISSB1 shows stress- and G3BP1-dependent nuclear efflux C_LIO_LISSB1 is upregulated and undergoes apical granulation in renal epithelial cells during renal ischemia-reperfusion injury C_LI

The Tar-DNA Binding Protein-43 C-terminal region, TDP43LC, has been previously shown to form amyloid-like fibrils with distinct folds in ALS and FTD. In both diseases, proteinaceous inclusions contain TDP43 C-terminal protein fragments as well as phosphorylated TDP43. Here, we use solution NMR to show that soluble phosphomimetic TDP43LC, P-TDP43LC, is structurally similar to wild-type TDP43LC. Disperse P-TDP43LC, like wild-type protein, contains a central helical region flanked by long disordered regions. Despite this similarity, our turbidity measurements, imaging, and kinetic assays show that P-TDP43LC has different aggregation behavior than wild-type protein. Using solid state NMR measurements we find that that phosphomimetic mutations alter the wild-type fibril conformation. Electrostatic repulsion from negatively charged sidechains, despite having little effect on the soluble proteins structure, perturbs amyloid-like fibril formation and selects for a different conformation in vitro. These results shed light on the structural role of TDP43LC phosphorylation in fibril formation in disease. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=104 SRC="FIGDIR/small/725298v1_ufig1.gif" ALT="Figure 1"> View larger version (16K): org.highwire.dtl.DTLVardef@1c63aforg.highwire.dtl.DTLVardef@1d48ed6org.highwire.dtl.DTLVardef@1ed8fd3org.highwire.dtl.DTLVardef@17d67a8_HPS_FORMAT_FIGEXP M_FIG C_FIG SynopsisPhosphomimetic mutations at ALS and FTD neurodegeneration-associated sites in an amyloid forming protein perturbs the aggregated structure compared to wild-type protein.

Bacterial genotoxins, Cytolethal Distending Toxin (CDT) and colibactin, cause DNA damage in intoxicated epithelial cells of the host. DNA damage influences {beta}-catenin signaling, altering its stability and nuclear translocation, potentially contributing to cancer development. Using non-transformed hepatocytes and cancer-derived intestinal and hepatic epithelial cell lines, we showed that CDT/CdtB induces the phosphorylation of {beta}-catenin at serine 552, along with the loss of {beta}-catenin from adherens junctions. This leads to the subsequent cytoplasmic accumulation and nuclear translocation of {beta}-catenin, ultimately driving TCF/LEF transcription, the crucial downstream event of Wnt/{beta}-catenin signaling, as well as the transcription of some {beta}-catenin target genes. Colibactin induces similar effects. MK-2206, a direct AKT inhibitor, and metformin, an AMP-activated protein kinase activator that indirectly inhibits AKT, both protected cells against various effects induced by CdtB exposure. These effects include {beta}-catenin phosphorylation at Ser552, the disassembly of cell-cell junctions and the subsequent nuclear accumulation of phosphorylated {beta}-catenin, leading to reduced TCF/LEF-mediated transcription. Additionally, MK-2206 and metformin protected from CdtB-induced epithelial-mesenchymal transition (EMT), including increased nuclear accumulation of SNAIL, enhanced matrix degradation and motility. Overall, these data show that infection with genotoxin-producing bacteria controls some EMT features through {beta}-catenin and AKT-dependent signaling.

The rapid expansion of human genomic data has revealed a large number of naturally occurring variants, creating a major challenge for functional annotation. The human metal transporter SLC39A8 (ZIP8) is a clinically important, promiscuous divalent metal transporter, yet most of its documented variants remain uncharacterized. Here, we developed a workflow to functionally evaluate ZIP8 variants by integrating laser ablation inductively coupled plasma time-of-flight mass spectrometry (LA-ICP-TOF-MS) with scaled-up cell-based transport assays. Using this method, we systematically analyzed 33 naturally occurring missense variants located in the extracellular domain (ECD) of ZIP8. The assay enables direct quantification of intracellular metal accumulation with substantially improved throughput ([~]150 samples per hour). Functional screening identified 14 potential pathogenic variants with significantly reduced transport activity. Comparison with computational predictions revealed a moderate correlation between activity and AlphaMissense pathogenicity scores (R2 = 0.423), while an error rate of [~]20% underscores the need for experimental validation. Flow cytometry analysis showed that most loss-of-function variants exhibit impaired trafficking of the protein to the cell surface possibly due to mutation-caused protein misfolding or instability. Structural mapping of activity-compromised variants, together with functional assessment of the ZIP8-ECD, highlights the importance of this domain in ZIP8 expression and intracellular trafficking. Together, this work establishes a scalable approach for functional screening of metal transporter variants and provides new insights into the structure-function relationships of ZIP8.

A viral exposure signature (VES) has been previously described predicting the development of Hepatocellular carcinoma (HCC) in at-risk patients. This has been achieved by a serological profiling of the viral infection history using a synthetic human virome including >100k epitopes (VirScan). In the present study we applied the same VirScan strategy to identify a differential serum binding pattern (DSBP) for classifying patients of different cancer types from healthy individuals. In particular, the healthy group included both age-matched (ADULTS) as well as elderly (ELDERS) individuals, the latter counting also nonagenarians and centenarians. The class comparison performed with serological data show DSBPs supporting class predictions, as confirmed by the receiver operating characteristic (ROC) curve analysis. Antibody responses supporting the class predictions are specific to peptides from persistent herpesviruses, acute-infecting viruses and, consistently in all comparisons, human respiratory syncytial virus (HRSV). Strikingly, the DSB of the ELDERS vs. CANCER comparison is characterized by higher titers in the healthy subjects; on the contrary, the DSB of the ADULTS vs. CANCER comparison is characterized by lower titers in the healthy subjects. Overall, the results show a differential serological binding pattern predicting healthy individuals (ADULTS or ELDERS) from patients with different types of cancer. Such results provide the first evidence suggesting a close link between anti-microbial immunity and cancer development. They may be of the highest relevance in terms of predictive, diagnostic and/or prognostic impact in oncology.

The functioning of proteins is intimately linked to the conformational states they sample within the native ensemble. Generating ensembles from a single static structure is therefore a research domain receiving considerable attention. In this application note, we introduce Hashi, a pipeline to rapidly generate realistic structural ensembles from the outputs of the structure-based Wako-Saito-Munoz Eaton (WSME) statistical mechanical model of protein folding. This approach relies on integrating the block WSME model outputs - strings of zeros and ones describing the conformational status of every residue over thousands or millions of microstates each assigned a statistical weight derived from physically grounded energy-entropy terms, and free energy profiles - with the RANCH module of the EOM (ensemble optimization method) from the ATSAS software suite, providing three-dimensional views of the structural ensembles within the model framework. It is applicable to a variety of single-chain monomeric systems with lengths ranging from 30 to 500 residues, including globular and repeat proteins. The generated structural ensembles can also be rank ordered according to their free energies within a given macrostate or a range of reaction coordinate values. Since the statistical weights of the WSME model microstates can be reweighted or calibrated with experiments, the ensembles shed light on not just the folding mechanism but also on the structural excursions that determine function and opening of otherwise buried binding pockets.

The insulin receptor samples multiple conformational states during ligand binding and activation, but the transient structural transitions connecting experimentally resolved receptor conformations remain poorly characterised. Here, we report a spontaneous opening event of insulin receptor site 1 observed during an unbiased molecular dynamics simulation initiated from the experimentally resolved singly insulin-bound IR1 receptor structure. The transition was characterised by separation of the L1 and FnIII-2 domains, rearrangement of site 1 contacts, and bending of the CT segment on the initially unoccupied receptor protomer. Structural comparison of the resulting conformation revealed similarity to the asymmetric IR2-A1 and IR2-A2 receptor states associated with hybrid insulin binding sites. Together, these findings suggest that conformations compatible with hybrid-site receptor states can emerge spontaneously from the intrinsic dynamics of the insulin receptor ectodomain, supporting a conformational selection model for receptor ligand engagement.

The human mitochondrial amidoxime reducing component 1 (mARC1) is a molybdenum-dependent enzyme whose protein-coding variants confer protection against common metabolic liver diseases. Whereas the frequent A165T variant acts largely through accelerated cellular degradation, the basis of protection by the rarer M187K variant remains obscure, as it has been reported that there are no differences between the "wild-type" and M187K variant protein in terms of cellular protein levels and localisation. Here, the crystal structure of the human mARC1 M187K variant, crystallised as a T4 lysozyme fusion after iterative micro-seeding, was determined at 1.63 [A] resolution. The variant structure is essentially superimposable with the previously reported "wild-type" and A165T structures, with pairwise root-mean-square deviations of 0.3-0.4 [A], and the pentacoordinated molybdenum cofactor is fully intact. Differential scanning fluorimetry across a broad pH range revealed only a modest, pH-dependent decrease in thermal stability associated with the exchange, most pronounced at alkaline pH. These data suggest that the disease-protective effect of M187K is unlikely to originate from gross structural rearrangement or active-site perturbation. SynopsisThe structure of the disease-relevant human mARC1 M187K variant was determined at 1.63 [A] resolution after iterative micro-seeding. The variant does not display relevant perturbations of the overall protein fold or active site structure, but differential scanning fluorimetry detects a pH-dependent decrease in overall protein stability associated with the M187K amino acid exchange.

Aptamers are single stranded DNA or RNA molecules selected for their high affinity and specificity to bind target molecules, similar to antibodies. They are commonly selected through the SELEX process, which involves the iterative exposure of a random sequence library to a target and retaining the sequences showing good binding properties. To improve Lyme disease detection, we propose designing aptamers that specifically bind to the CspZ protein on the surface of Borrelia burgdorferi, the bacterium responsible for the disease. Starting with a SELEX process consisting of thirteen rounds, from which selected in vitro sequence candidates have emerged, we aim to propose a holistic process that selects in silico new sequence candidates that are further validated experimentally. Our approach relies on 1) using Machine Learning (ML) techniques, specifically a Restricted Boltzmann Machine (RBM), to digitally replicate the last round of the SELEX process, 2) integrating insights from text analysis methods, such as word2vec and n-grams, into the RBM model trained on the final-round SELEX dataset to represent and compare newly generated sequences with in vitro candidates, 3) selecting in silico sequences with strong potential to bind to CspZ protein, 4) experimentally validating the selected in silico sequences of step 3. Our holistic approach combines biological insights with statistical models to improve the efficiency and outcome of the SELEX process. We enhance the RBM model, designed to replicate the distribution of the final SELEX round, by integrating geometric representations of sequences, which is especially advantageous when dealing with limited datasets relative to the vast sequence space. In addition, it provides in silico sequence candidates with strong binding properties.

Despite progress in predicting protein structures, how proteins arrive at their native state remains a subject of continuous debate. We present a single molecule force spectroscopy study of the unfolding and refolding intermediates of the conserved, diverse, and ancient Rossmann2x2 fold ({beta}12{beta}34{beta}56{beta}78). By inserting glycines at different locations in the protein, we can follow in real time and annotate its unfolding and refolding intermediates. This protein folds along a single reversible pathway involving the ordered and sequential organization of discrete and cooperative folding units or foldons: unfolded {rightleftarrows} {beta}12{beta}3 {rightleftarrows} {beta}12{beta}34{beta}5 {rightleftarrows} {beta}12{beta}34{beta}56{beta}7 {rightleftarrows} {beta}12{beta}34{beta}56{beta}78. This strict order results from the formation of an autonomously folding unit (primary foldon) and the subsequent organization of elements (secondary foldons) whose stability depends on their interactions with previously organized ones.

The molten globule (MG) state is an intermediate in the unfolding pathway of proteins, typically triggered by denaturing agents such as urea, extreme pH, high pressure, or heat. The microscopic details of such states are far from understood. Here we study the MG states in protein Hen Egg-White Lysozyme (PDB ID: 1AKI) using microscopic constant pH molecular dynamics (CpHMD) simulations and experiments across a wide pH range. We observe that the titratable residues act as key drivers of conformational fluctuations, promoting the emergence of MG states at extreme pH. These states display partial unfolding, and small global structural changes (< 7% deviation). Hydration around the fluctuating acidic residues shows reduced water density and weakened hydrogen bonding at low pH. At high pH, hydration around acidic residues increases relative to pH = 7, whereas hydration around basic residues decreases. The translational and rotational dynamics of hydration water also exhibit pronounced pH dependence: the translational diffusion coefficient (Dtrans) increases linearly with decrease in pH in acidic medium and increases linearly with increasing pH in the basic regime. The rotational diffusion (Drot) shows similar dependencies on pH except a break at pH {approx} 4 corresponding to acidic residue pKa values. Our results may be useful to identify ligand binding of lysozyme in extreme pH conditions.

CD239, also known as Lutheran blood group glycoprotein (Lu) or basal cell adhesion molecule (B-CAM), is a transmembrane protein belonging to the immunoglobulin superfamily (IgSF). CD239 serves as a specific receptor for laminin 5, a subunit of laminin-511/-521, which are major components of renal basement membranes. A previous study of another group demonstrated that CD239-null mice are healthy and develop normally. Although no alteration in renal function was observed, most glomeruli in the mutant kidneys exhibited morphological abnormalities. In this study, we investigated the role of CD239 in renal tubules. We examined the distribution of CD239 using renal tubule-specific markers. CD239 was localized to the Henles loop, distal tubule, and collecting duct, but not to the proximal tubule. Next, we analyzed the localization of renal tubular molecules in CD239-null mice. The localization of uromodulin (UMOD) and Na-K-Cl cotransporter (NKCC2) was disrupted in the distal tubules lacking CD239, suggesting that CD239 plays a role in maintaining the polarity of renal epithelial cells. Furthermore, to examine the stability of the distal tubules, CD239-null mice were subjected to chronic kidney disease (CKD) using an adenine-rich diet. Blood analysis revealed that CD239-null mice fed an adenine diet readily developed CKD. Adenine-fed null mice exhibited more marked histological injury along the distal tubules compared to that by controls. These results indicate that CD239 is essential not only for maintaining cellular polarity but also for ensuring the stability of the distal tubules. Although CD239-null humans exhibit no obvious associated pathology, it could be a predisposition to CKD.

BackgroundHomozygous deletion of the methylthioadenosine phosphorylase (MTAP) gene is a frequent genetic alteration in cancer. MTAP, which creates adenine from 5-methylthioadenosine (MTA), is constitutively expressed in all tissues throughout the body. Previously, we described a novel strategy to specifically target MTAP-deleted cancer cells by combining the antipurine prodrug 2-fluoroadenine (2FA) with MTA. In vitro, this combination efficiently killed MTAP- cancer cells, but in vivo the combination was much less effective in vivo. Here, we explored the role of xanthine oxidase (XO) in this process. Materials and MethodsVarious combinations of 2FA, MTA, and the xanthine oxidase inhibitor febuxostat (FX) were tested in various cancer cell lines grown in vitro and in mice. LC-MS/MS was used to examine the levels and ratio of intracellular 2-FA-containing nucleotides compared to adenine-containing nucleotides. Results and conclusionsThe treatment of cells with 2FA+MTA in vitro resulted in much higher 2FANP/ANP ratios than the same treatment in vivo. The addition of XO to culture media in vitro effectively abolished the killing by 2FA, and this effect was fully reversed by the addition of febuxostat (FX), a xanthine oxidase inhibitor. In vivo, the addition of FX to 2FA results in increased cell killing and toxicity and a 1000% increase in the amount of 2FA converted to 2-FA-monophosphate (2FAMP). Xenograft studies using MTAP- HT1080 and MiaPaCa-2 cell lines have shown that a 2FA/MTA/FX cocktail can cause tumor regression in vivo. These studies suggest that the combination of 2FA/MTA/FX should be explored as a treatment for MTAP- cancer.

HIV-1 buds from infected cells as immature virion particles with a scattered envelope glycoprotein (Env) distribution on their envelope. It then undergoes maturation, during which the viral protease cleaves the Gag polyprotein at multiple sites, leading to structural reorganization of the viral particle and lateral redistribution of Env proteins, ultimately rendering the virion infectious. However, the underlying mechanism of maturation-induced Env reorganization remains elusive. In this study, we combine microsecond-long all-atom (AA), bottom-up coarse-grained (CG) molecular dynamics simulations, and diffusion model-based backmapping to investigate the structural organization and key interactions of Env in viral membranes. AA simulations of fully glycosylated Env embedded in HIV-1 mimetic asymmetric bilayers were first performed to characterize its conformational dynamics and Env-lipid interactions. We then developed a bottom-up CG model of glycosylated Env from that AA data and simulated the mature HIV-1 virion envelope containing multiple Env proteins. The CG simulations predict that Env proteins form clusters through interactions mediated by the cytoplasmic tail domain (CTD) and adopt diverse tilted conformations within these clusters. These CG simulations were then backmapped to AA resolution and further AA simulations were carried out to identify, in detail, the specific interacting residues in the Env clusters. Additionally, analysis of epitope accessibility shows that broadly neutralizing antibodies (bnAbs) targeting the V1/V2 and V3 loops may efficiently interact with Env clusters on the mature virion surface. Together, these results provide a molecular mechanism for Env oligomerization during viral maturation and offer new insights into the accessibility of bnAb epitopes on Env clusters.

We present Pathway Representation via Intrinsic Structural Medoids (PRISM), a state-aware framework for clustering pathways from molecular dynamics simulations of biomolecular transitions. In PRISM, each pathway is mapped to a small set of structural medoids obtained via a deterministic k-means clustering scheme. Pairwise pathway dissimilarities are computed using a weighted average Hausdorff distance between these representative sets, effectively capturing mean nearest-neighbor structural deviations while reducing sensitivity to outliers. Hierarchical agglomerative clustering of the resulting dissimilarity matrix defines pathway families. We evaluate PRISM across three biomolecular transitions of increasing complexity: alanine dipeptide C7eq [-&gt;] C7ax isomerization, adenylate kinase opening, and HIF-2 PAS-B ligand unbinding. PRISM consistently yields robust cluster assignments, with medoids faithfully representing distinct conformational states. By combining a state-based description with robust geometric dissimilarities, PRISM provides a scalable framework for organizing complex transition pathways.